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Glass is a microscopically disordered, solid form of matter that results when a fluid is cooled or 
compressed in such a fashion that it does not crystallise. Almost all types of materials are capable 
of glass formation polymers, metal alloys, and molten salts, to name a few. Given such diversity, 
organising principles which systematise data concerning glass formation are invaluable. One such 
principle is the classification of glass formers according to their fragility Fragility measures the 
rapidity with which a liquid's properties such as viscosity change as the glassy state is approached. 
Although the relationship between features of the energy landscape of a glass former, its configurational 
entropy and fragility have been analysed previously (e. g., [^), an understanding of the origins 
of fragility in these features is far from being well established. Results for a model liquid, whose 
fragility depends on its bulk density, are presented in this letter. Analysis of the relationship between 
fragility and quantitative measures of the energy landscape (the complicated dependence of energy on 
configuration) reveal that the fragility depends on changes in the vibrational properties of individual 
energy basins, in addition to the total number of such basins present, and their spread in energy. 
A thermodynamic expression for fragility is derived, which is in quantitative agreement with kinetic 
fragilities obtained from the liquid's diffusivity. 

Glass forming liquids grow increasingly viscous upon cooling, till the viscosity becomes so large that they fail to 
flow on experimental time scales. While remaining microscopically disordered like a liquid, they manifest mechanical 
properties of a solid; i. e., they transform to a glass. By convention, the glass transition temperature Tg is that where 
the viscosity reaches a value of 10^^ Pa sec. The approach to this large viscosity, however, differs from one liquid 
to another. When displayed in an Arrhenius plot [log(viscosity) vs. inverse temperature l/T], some liquids (e. g., 
silica) show a steady, linear increase, while others display a much steeper dependence on l/T |^,^, as illustrated in 
Fig. 1(a) (inset). The former are strong liquids, and the latter, fragile. This range of behaviours is implicit in the 
Vogel-Fulcher-Tammann-Hesse (VFT) form, observed to describe the T dependence of viscosity (as well as diffusivity 
and relaxation times) in many glass formers. The VFT relation may be written as 
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where Tq is the temperature of apparent divergence of viscosity, and Ky ft is a material specific parameter quantifying 
the kinetic fragility; more fragile liquids have larger KypT values. This behaviour of viscosity correlates with the jump 
in heat capacity at the glass transition - more fragile the liquid, sharper and bigger the jump. Rationalisation of this 
correlation comes from the Adam-Gibbs relationship Q which predicts a dependence of viscosity on the configurational 
entropy Sc (described later) of the liquid: 
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which results in the VFT relation if Sc has the form TSc = Kag{T/Tk — 1) (2^, where Tk is the ideal glass transition 
or Kauzmann temperature below which Sc ~ 0. The jump in heat capacity is proportional to Kaq., and is larger 
for more fragile glass formers. 

Fragility is analysed here for a binary mixture of particles interacting via the Lennard- Jones potential, widely studied 
as a model glass former [|]-^. The T dependence of diffusivity D (to characterise dynamics) and configurational 
entropy are obtained for a range of bulk densities p. The diffusivities D{p,T), are shown in Fig. 1(a) in a 'fragility 
plot': — log[-D(p, T)] vs. Tr/T, where Tr is a reference temperature (akin to Tg in the usual fragility plot) at which D 
at each density reaches a fixed, small value. By comparison with the inset, the liquid is seen to become more fragile as 
its density increases. Fig. 1(b) shows the kinetic fragility index, Kvft, obtained from VFT (Eq.(|^)) fits to diffusivity 
data at each density. The density dependence of Ky pT quantifies the increase of fragility with density, apparent from 
data in Fig. 1(a). 

A thermodynamic explanation of fragility, as discussed above (also |l^,|ll[]), lies in explaining the rapidity of change 
with T of the configurational entropy, Sc- In the inherent structure (IS) formalism |12jl^ employed here, the con- 



figuration space of a liquid is divided into basins of local potential energy minima ('inherent structures'). Sc at a 
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given T derives from the multiplicity of local potential energy minima sampled by the liquid at that T, and equals 
the difference of the total entropy and the vibrational entropy of typical basins sampled, Smb'. 

Sc{p,T) = Stotai{p,T) - Syib{p,T). (3) 

For all the results presented here, the basin entropy Syib is calculated by approximating each basin as a harmonic 
well (valid at sufficiently low T) 0-^,0. In this approximation, the entropy of a given basin, arising from vibrational 
motion within the basin (the suffix vib emphasizes this point) is given by 

37V , 

where Vi are the vibrational frequencies of the given basin, ks is Boltzmann's constant and h is Plank's constant. 
The basin free energy is Fyib — Uyib — TSyib, where Uyib — ZNksT is the internal energy. Apart from the explicit 
temperature dependence of Syib above, an implicit dependence also exists because vibrational frequencies change 
from one basin to another, and, different basins are sampled at different temperatures. From the expression above, 
it is clear that the entropy difference between any two basins arises solely from the difference in their vibrational 
frequencies, and is independent of temperature. 

The configurational entropy is evaluated in this letter using two methods, referred to as the thermodynamic integra- 
tion (TI) method, and the potential energy landscape (PEL) method, respectively. In both cases, inherent structures 
are obtained by performing local potential energy minimisation for a subset of configurations sampled by a liquid 
at each p and T . Vibrational frequencies are calculated for each inherent structure, and the basin entropy Syib is 
obtained from Eq. (^). In the TI method (described in detail elsewhere see also caption of Fig. 2), the total 

entropy S total is obtained using pressure and internal energy values from simulations. Eq. (|^) is used to obtain Sc- In 
contrast, the PEL method is based on constructing the distribution of inherent structures in order to calculate Stotai 
and Sc, as described later. 

Fig. 2(a) shows the Adam-Gibbs plot, where logD is plotted against [TSc)^^. The Adam-Gibbs relation (Eq. (||)) 
is seen to apply well at each density, as also observed previously for the hard sphere liquid jl^ and water Fig. 
2(b) shows (as lines) the T dependence of TSc (from TI). All curves are nearly linear, but with slight increases of slope 
at lower T . The slopes of TSc (hi the T range where diffusivities are measured), K^q, are used as the quantitative 
index of thermodynamic fragility, and are plotted in Fig. 1(b). are nearly the same as the kinetic fragility index 

KvFT, thereby validating the idea that the rate of increase of T Sc does indeed determine the fragility. 

However, this analysis provides no insight into how the fragility relates to features of the energy landscape. To 
this end, the distribution of energy minima, and properties of individual basins, are now examined, at fixed density 
p. The probability distribution P{^, T) that an inherent structure of energy $ is sampled at temperature T depends 
on (a) the energy (b) the structure of the basin, contained in the basin free energy Fyib{^,T) and (c) the number 
density of inherent structures with energy [|2|jl|J^,|l|,|l^ : 

P($, T) = r!($) exp [-/?($ + Fy.bi^, T))] /Qn{p, T), (5) 

where Qn, the partition function, normalises P and is given by Qn — J P{^,T) over all possible <&. The configu- 
rational entropy density is defined by 5c($) = fc^ lnr2($) and is related to the T dependent configurational entropy by 
Sc{T) = / d$ 5c($)P(5>,r). The quantities P($,T), Fyib and Qn are obtained from simulations |0,^|j|Jl^, and are 
used to invert the relation in Eq. ^ to obtain the configurational entropy density iSc($). The same 5c($) estimates 
must result (though not for the same range of <&) for different T, when the harmonic assumption is valid, which is 
indeed seen to be the case from data shown in Fig. 3(a). With increasing density, the maximum 5c($) value reached 
decreases, and the distribution becomes broader. It has been argued ||l^,|l^ that the distribution ^{^) should be 
Gaussian (equivalently, Sd^) an inverted parabola). Data in Fig. 3(a) are fully consistent with this expectation, and 
are fitted to the form, 

Scm ($-^o)^ ... 

Nkb (J 

where a is the height of the parabola and determines the total number of configurational states, i. e. energy minima 
(the total number is proportional to exp{aN) ), $0 and respectively define the mean and the variance of the 
distribution. The fits are displayed in Fig. 3(a) and the fit parameters listed in the figure caption. Speedy |^ has 
considered a model of the thermodynamics of a glass former with a Gaussian distribution of configurational states, 
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and the basin specific heat that does not depend on the internal parameter (here the IS energy). He predicts that 
the fragiUty of a liquid should be proportional to the logarithm of the number of states, i. e. on a. From Fig. 3(a) 
it is clear that a values, and hence the number of states, decrease with increasing density p. But the fragilities (Fig. 
1(b)) show the opposite trend. However, fragility has an additional dependence (implicit in ref. [^) on the variance 
of the distribution a (see equation (6)). 

Further, the variation of vibration frequencies, and hence the basin entropy, from basin to basin, also contributes 
to the fragility of the liquid. Fig. 3(b) displays the T independent (see above) basin entropy difference at a given 
energy $ with respect to the value at <I>o, ^Smb{^) = Svib{^,T) — Syib{^Q,T). At all densities, ASmb{^) displays 
a clear $ dependence, which is inconsistent with assumptions in ||^ . Although such basin dependence of vibrational 
frequencies and hence the basin entropy have been discussed previously |^-^ , with reference to both experimental 
data and model calculations, its implications on a liquid's fragility have not clearly been established. Fig. 3(b) 
shows that basins at higher energy (which are occupied at higher temperature) have lower basin entropies, or, higher 
frequencies (see Eq. (^)). This trend, contrary to what one expects for systems studied experimentally, arises because 
the constant density model system here cannot expand when T is raised. 

The dependence of ASvib{^) on <I> is to a very good extent linear. Hence, for each density, ASvib is fit to the form 
ASvibi^) = SS — $o) with fit parameters listed in the caption of Fig. 3. The basin free energy, which can be 
written as Fyib{^,T) = Fyib{^o,T) — TSS{^ — $o)j is also linear in $. With a Gaussian form for Q, the partition 
function Qn = J d^P{^,T) [P as in Eq. (|^)) is easily evaluated, and results in the following predictions for the T 
dependence of the average IS energy, and the configurational entropy: 

where ^l^^ = *o + and 
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where Tk — <j{2NkB\/a + cr5S)~^ is the ideal glass transition temperature where Sc{Tk) = 0. 

The form above for Sc results in the VFT relation, if K^q^{T) is constant. This is not quite the case, but the 
1/T term becomes rapidly irrelevant for T > Tk and K^q^{T) approaches a constant. Further, TSc values plotted 
(as points) in Fig. 2(b) are in extremely good agreement with the ones from thermodynamic integration, including 
deviations from linearity for small T /Tk- A check of the overall consistency of the model described here is made 
in Fig. 3(c), which shows a scaled plot of < $ > (T) — '^q''^''^ vs. 2NkBT • From Eq. (^, one expects $(r) data 
for all densities to collapse onto a straight line of (negative) unit slope, with expected deviations at high T due to 
basin anharmonicity |p^ . Data in Fig. 3(c) clearly meet this expectation, although with noticeable deviations for 
p — 1.35po, as also in Fig. 2(b). 

The expression for K^q^ in Eq. (^, the thermodynamic fragility index obtained from energy landscape parameters, 
shows that fragility depends not only on the multiplicity of states (a) but also on their spread (ct) , and how much the 
basin entropy changes from the lowest to the highest energy basins sampled by the liquid. The extent of this change is 
quantified by parameters a and 6S (listed in the caption of Fig. 3). Fig. 1(b) shows the average slopes of TSc (from 
Eq. (^) above, in the T range where diffusivities are measured), K^q^ , which are nearly the same as Kj^^, and in 
very good agreement with the kinetic fragility index Kvft- This agreement establishes the quantitative relationship 
between fragility and the energy landscape of the liquid since K^q^ is obtained solely in terms of a quantitative 
description of the liquid's energy landscape. 
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Figure 1: Fragility Plot of difFusivitics, kinetic and thermodynamic fragility indices. Results shown are obtained 
from molecular dynamics simulations of 204 A type and 52 B type particles interacting via the Lennard-Jones (LJ) 
potential. Argon units are used for A type particles: LJ parameter eaa = ks ^ 119.%K, uaa — 0.3405 x 10~^m, 
VTLA = 6.6337 X 10~^^fc(7, and e^s/eAA — 1-5, ^bb/^aa = 0.5, ctab/caa = 0.8, and OBBl^yAA — 0.88, mB/mA — 1. 
Densities are reported in units of po = 2.53 x lO^^m"'^ or 1.678gcm~^ . The liquid-gas critical point is located roughly 
at Tc = 130i4', pc = 0.416po- (a) Diffusivities D from molecular dynamics simulations (details of the simulations 
are as in |^) are displayed in a 'fragility plot' for 5 densities. The reference temperature Tr is chosen such that 
D{Tr) — 2.44 X 10~*cm^/s, which is slightly below the lowest D values measured. The VFT extrapolation is used 
to locate Tr- Fit values Do (preexponent) are in the range of 0.87 to 2.0 x 10~'*cm^/s, close to the experimentally 
observed values, (b) Fragility index K obtained (i) from VFT fits to diffusivities D, and from the configurational 
entropy obtained (ii) from thermodynamic integration ('TI'), and (iii) from analysis of the potential energy landscape 
('PEL'). Dimensionless thermodynamic fragility indices are obtained by dividing if Jq, K^q'" by NeAA- Comparison 
with data in [0 shows that at the low density end, the fragility of the model liquid compares with that of considerably 
fragile liquids such as toluene, orthoterphenyl and salol, while at the high density end, the fragility is extremely high. 
Temperatures of vanishing diffusivity To from VFT fits are 18.76i4:, 24.98is:, 35.97if , 47.86if and 75.75/f for 
p/ po = 1.1, 1.15, 1.2, 1.25 and 1.35 respectively. 

Figure 2: Adam-Gibbs plot, configurational entropy vs. temperature, (a) The 'Adam-Gibbs' plot displaying 
log(D) vs. [T S'c)~^. T Sc is expressed in units of caa- The Adam-Gibbs expectation of a linear dependence is 
well satisfied at each density, (b) Configurational entropy obtained from thermodynamic integration ('TI') (lines, not 
labeled), and from analysis of the potential energy landscape (points, labeled as 'PEL'), plotted against temperature 
T scaled to the Kauzmann temperature Tr- obtained from thermodynamic integration. In the TI method, the total 
entropy Stotai of the hquid is obtained from its total free energy Ftotai, from Ftotai = U total - TStotaU where Utotai is 
the internal energy. Thermodynamic identity P — p^ ( ^^*°dp' ) where P is the pressure, is used to integrate P from 

simulations to obtain Ftotai with respect to the known ideal gas free energy. The identity Utotai — ( ^^biJ't'^ ) used 

to integrate at fixed density Utotai from simulations to obtain Ftotai at any desired T. Extrapolations to low T beyond 
the range of simulations utilize the form for the potential energies E, E = Eo + EiT^^, where E2 ~ 3/5 (see [0-|)). 
Basin entropies S^b £^re obtained by averaging over typical inherent structures at each T, whose vibrational frequencies 
are obtained numerically by diagonalising the Hessian (matrix of second derivatives of the potential energy). The 
plot covers at each density the range of T where simulations have been performed. Scaling temperatures , the 
Kauzmann or ideal glass transition temperature where Sc = 0, are: 2Q.AAK, 26.83-fir, 34.35^, 'i2.49K, and 62.A7K, 
for p/pa = 1.1, 1.15, 1.2, 1.25 and 1.35 respectively. T^^^ obtained from the PEL method are, in the same order, 
20.24X, 26.18X, M.7AK, and 62.35X. 

Figure 3: Distribution of energy minima, basin entropy, and temperature dependence of inherent structure energies. 

(a) The configurational entropy density 5c($) plotted against inherent structure energy $ minus the 'Kauzmann' 
energy $mm below which 5c ('I') — 0. is the lower root of 771^ = a — ('^~'|'°) — q Pqj- each density, 
data from 9 temperatures (on average) are shown (smaller symbols) which collapse onto each other, validating the 
harmonic approximation for calculating basin entropies (see text). Also shown are the averages of these curves 
(larger symbols), for whose calculation a threshold of 20 samples (out of a total of 1000) is imposed in each bin. 

(b) The T-independent component of the basin entropy AS'^if,(<I>) vs. $ — $0, where $0 is the value where 5c(<I') is 
maximum. ASyib{^) curves are approximately linear. The T independence of AS'„if,($) is an artefact of the classical 
mechanical treatment employed here, but a quantum mechanical calculation of the basin entropy does not alter 
the picture fundamentally (unpublished results), (c) Scaled plot of the average IS energies sampled in simulations, 
which are expected to collapse on to the straight line shown at low temperatures (right side of the graph), if the 
assumption of a Gaussian configurational entropy density Sc{^) is valid (see text). The observed data collapse 
validates the Gaussian form used for Sc{^). Deviations at high temperatures are due to the breakdown of harmonicity 
of basins. Fit parameters a, (j)o/N€AA, cr/NeAA and 6SeAA/kB are, {p — l.lpo): 1.010, — 6.586, 0.303, — 0.721, 
(p = 1.15po): 0.957, - 6.684, 0.374, - 0.652, (p = l.2po): 0.933, - 6.683, 0.489, - 0.500, [p = 1.25po): 
0.927, - 6.564,0.642, - 0.314, and (p = 1.35po): 0.849, - 6.088,0.873, - 0.188. 
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